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1 Introduction 



Over the years, the exact solution of A^-body problems has attracted consid- 
erable attention because of its possible relevance in statistical mechanics as 
well as in atomic, nuclear and gravitational many-body problems. Whereas 
exact solution of several iV-body problems in one dimension is known by now 
[1,2], a class of exact solutions including the bosonic ground state have been 
obtained for several A^-body problems in higher dimensions when they are 
interacting via a harmonic oscillator potential [3,4,5,6,7,8]. Clearly it is of 
considerable interest to discover other exactly solvable A^-body problems in 
one as well as in higher dimensions. 

The purpose of this paper is to show that a class of exact solutions in- 
cluding the bosonic ground state of all these A'"-body problems in higher 
dimensions can also be obtained in case they are interacting via an A^-body 
potential of the form 

V(ri,r2,...,riv) = ^= (1) 

In this context, 1 may add that recently 1 have already obtained a class of 
exact solutions of the A^-anyon problem (in two dimensions) in case they are 
interacting via this A^-body potential [9] . Further, last year, I also obtained 
the complete bound state spectrum of the A^-particle problem in one dimen- 
sion in case they are interacting via a variant of the above potential as given 
by [10,11,12] 

V{xx,X2,...,xn)^ — =====. (2) 

Subsequently, Gurappa et al. [13] showed that the complete bound state 
spectrum can also be obtained in case an A^-body potential of the form 
P'^/J2i<j{^i — XjY is added either to the oscillator or to the A'"-body potential 
(2). 

Based on all these results, I conjecture that whenever an A'^-body problem 
is solvable in case these A^ particles are interacting via an external one body 
( or pairwise) oscillator potential, the same A^-body problem is also solvable 
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in case they interact via the A^-body potential as given by eq. (|I]) or by 



I further conjecture that in either case, one can also add an iV-body potential 
of the form 

r(ri,r2, ...,r^) = -— — -, (4) 



or its variant 



y(ri,r2,...,riv) = -— — 2 (5) 

and the problem is still solvable except that the degeneracy in the spectrum 
is now much reduced. Clearly, it is necessary to examine the known solvable 
A^-body problems with the oscillator potential and check if these conjectures 
are valid or not. This is precisely what I propose to do in this paper. 

The plan of the paper is as follows. In Sec. II, I discuss the A^-body 
Calogero-Marchioro [3] model in D-dimensions [4] and show that a la the 
oscillator case, even for the A^-body potential as given by eq. (^, one can 
obtain some exact eigenstates including the (bosonic) ground state. I also 
consider the Sutherland variant of the problem and as in the oscillator case, 
I obtain the exact (bosonic) ground state when the A^-bodies are interacting 
via the A^-body potential as given by eq. (1). In Sec. Ill, I discuss the two- 
dimensional model of Murthy et al. [5] in which they have obtained some 
exact eigenstates including the (bosonic) ground state all of which show novel 
correlations. I show that similar exact solutions with novel correlations can 
also be obtained in case one replaces the one-body oscillator potential by 
the potential as given in eq. (1). I also discuss the two-body problem in 
great detail and show that as in [6] , in this case too it is completely solvable. 
In Sec. IV, I consider the D-dimensional generalization [7] of Murthy et al's 
model [5,6] and show that exact solutions including the (bosonic) ground 
state can again be obtained in case they are interacting via the A^-body 
potential (1). In Sec.V, I consider a Calogero type model in D-dimensions 
[8] which has only two-body interactions and show that one can obtain some 
exact eigenstates including the (bosonic) ground state in case the A^-bodies 
are interacting via the A^-body potential (1). Finally, in Sec. VI, I obtain the 
entire discrete bosonic spectrum in case the A^-particles are interacting in 
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/^-dimensions purely via the A^-body potential as given by (^. Further, in 
all the above cases I show that the discrete spectrum can still be obtained 
even if one adds the iV-body potential as given by eq. (^) (or (|)) to the 
oscillator or our A^-body potential except that the degeneracy in the discrete 
spectrum is now much reduced. 



2 Calogero-Marchioro Model With A^-body 
Potential 

Long time ago, in an effort to generalize the original Calogero model [2] 
to dimensions higher than one, Calogero and Marchioro [3] considered a 
model in three dimensions in which the N-particles are interacting via the 
two-body and the three-body inverse square interactions as well as by the 
pairwise harmonic oscillator potential. They were able to obtain some exact 
eigenstates including the bosonic ground state of the system. Recently, we 
[4] have been able to generalize their work to arbitrary space dimensions. In 
particular we considered the following A^-body Hamiltonian in D-dimensions 

rT_ VY72 . gh^^ I Gh^^ Vki-Vkj rnu^^ , 

and showed that some eigenstates including the bosonic ground state for the 
system are given by 



with energy 



2 V 27V ^ ^ ' v 27V 



En. = J Y(2n, + Tz) + 1) , = 0, 1, 2, ... . (8) 



Here is the D-dimensional position vector of the z'th particle and rij = 
Fj — Tj denotes the relative separation of the i'th and j'th particles while r^j 
denotes its magnitude. In writing this exact solution, we have scaled all dis- 
tances Tj \J ^'^d the energy is measured in units of Tluj. Throughout 
this paper, whenever we discuss the exact solutions for the oscillator poten- 
tial, we shall always be working with these scaled coordinates and the energy 
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will always be measured in units of huj. The parameters g and G character- 
ize the (dimension-less) two-body and three-body coupling constants respec- 
tively while denotes the Laguerre polynomial. Here g > —1/2 to stop 
the fall to the origin. Further, Ad and Td are two parameters determined in 
terms of the parameters of the Hamiltonian by 



Ad = VG 



[D - 2)2 + - (/} - 2) 



r 



D 



D{N - I) - 2 + AdN{N - I] 



(9) 



(10) 



Let us now consider the same many body problem as given by eq. (|^) 
but with the oscillator potential being replaced by the A^-body potential as 
given by eq. (^). Throughout this paper, whenever we talk of the exact 
solutions with the A^-body potential as given by eqs. (|^) or (^), we shall 
rescale all distances h^ri/me^ and measure energy in units of me^/h^ 
so that m, e, h are all scaled away. 

On substituting the ansatz 



^2 \Ad/2 



(11) 



in the Schrodinger equation for the potential (H) one obtains 



b"{p) + {2Tn + l)0'(p) + ( 



N 



where F^) is as given by eq. ( [TI1| ) while 

1 



P 



On further substituting 



2p I E |)0(p) = 



(f){p) = exp{-^2 I E \p)x{p) 
it is easily shown that x(p) satisfies the equation 



(12) 



(13) 



(14) 



PX"(P) + (2F^+1-2J2 I E \p)x'{p) + [^~{2Tn+l)j2 \ E \]x{p) = (15) 



N 
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whose solution is a Laguerre polynomial 



x{p) = I^l-{2^2\E\p) (16) 
with the corresponding energy eigenvalues being (in units of me^/h'^) 

= -2iV(n, + rB + i)2 • 
Several comments are in order at this stage. 

1. For D = 1 this expression agrees with the one derived by us earlier 
[10]. 

2. For rir = 0, the eigenfunction 

^0 = (n 4)^"^' I Eo \p) (18) 

i<j 

has no nodes besides those coming from the singular centrifugal poten- 
tial which forces the eigenfunction to vanish whenever the coordinates 
of any two particles coincide. Thus ipo corresponds to the ground state 
of the system with the corresponding ground state energy being as 
given by eq. ([TTD with = 0. For > 0, we have radial excitations 
over the ground state. 

3. As (7 ^ 0, we see from eq. (|^) that also G ^ and the wave func- 
tion (|18D becomes the ground state eigenfunction of the hyperspherical 
"Coulomb problem" in D-dimensions without centrifugal barrier and 
with Bose statistics. Thus the situation is different from the one di- 
mensional problem [10] where as g 0, the eigenfunction is the ground 
state of the "Coulomb problem" but with Fermi statistics. We shall 
in fact see that in all the higher dimensional many body problems 
{D > 1), unlike the one dimensional the coupling is switched 
off, the eigenfunction corresponds to that of Bose statistics. It is not 
clear whether this difference in statistics between the one and the higher 
dimensions has any deeper physical significance. 

4. The A^-body problem is still solvable if apart from our A^-body potential 
d^) we also add the potential (^ to the Hamiltonian (^. In this case, 
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on substituting the ansatz (|TT]) in the Schrodinger equation for the 
combined potential yields 



^"W + ^^^^'(P) + - - 2 I £ l)^(p) =0. (19) 

On following the steps as given by eqs. (|T3]) to ( ]T7| ) it is easily shown 
that the exact eigenstates are given by (in units of me^/fi^) 



0(p) = p(^-r-)exp(-V^2 I E \p)Llli2^2 \ E \p) (21) 

where 



^ = ^Tl + 6yN. (22) 

5. In the same way, the A^-body problem with the oscillator potential as 
given by eq. (^ is also solvable in case we also add the A^-body potential 
(^) to it. In particular, it is easily shown that the corresponding exact 
eigenstates are 

^n. = J^[2n, + l + 7] (23) 



i<j ' i<j ' i<j 

It may be added here that following Calogero and Marchioro [3], we can 
also obtain a subset of the (spin-less) fermionic eigenf unctions in the special 
case of = 3 and 4 when these particles are interacting via the A^-body 
potential as given by eq. (^. In particular, it is easily shown that for A^ = 3, 
a set of completely antisymmetric eigenstates in three space dimensions are 



i<j 



^n. = (ri2 X r23)(n4)^^'e^^^(-V2 I E \p)L'/A2J2 \ E \p) (25) 
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where 

_ 9. r I 1 

F = 3/ + 3 . (27) 




Generalization to D dimensions is straight forward. 

On the other hand for = 4, a set of completely anti-symmetrical eigen- 
states in three space dimensions are given by 



i^n. = [(ri2 X r23) ■ r,,]{l[viy/'expi-p \ E \p)L'Z' {2p \ E \p) (28) 



8(n, + F' + i^2 (29) 



2' 



where f is as given by eq. (pT]) while F' = 6f + 5. Again, generalization to 
D dimensions is straight forward. 

I might add here that a la the oscillator case [4], we can also obtain the 
exact ground state of the corresponding Sutherland variant of the problem. 
In particular, consider the Hamiltonian as given by eq. (^ but with the 
oscillator potential being replaced by the potential as given by eq. (|l|). It is 
easily shown that the exact ground state of the system is given by 



^0 = (11 4)^"^' exp{-j2 I E I Y: r?) (30) 

where 



i<j 



iN-ini + NAnY' ^^^^ 
One of the unsolved problem is whether one can map this problem (at least 
in some specific dimension D) to some random matrix problem and obtain 
exact results for the corresponding many-body theory. In this context it may 
be noted that for g = 2, the corresponding oscillator problem (see eq. (H)) 
in two dimensions has been shown to be connected to the random matrix 
problem for the complex matrices [4]. 



3 Novel Correlations With an A^-body Poten- 
tial 

In a recent paper, Murthy et al. [5] have proposed a model in two dimensions 
with two-body and three-body interactions. They were able to obtain the 
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bosonic ground state and a class of excited states of the system exactly by 
adding an external (one body) harmonic oscillator potential. The model was 
constructed in such a way that the solutions have a novel correlation of the 
from 

Xij = Xit/j Xjt/i (32) 

built into them. Note that unlike the Jastrow-Laughlin form, Xij is a pseudo- 
scalar. However, unlike the Laughlin type of correlation, this correlation is 
not translationally invariant unless the radial degrees of freedom are frozen. I 
now show that the bosonic ground state and the radial excitations over it can 
also be obtained in case the oscillator potential is replaced by the A^-body 
potential as given by eq. (0). Further, I also obtain the complete solution of 
the two-body problem. Notice that the two-body problem is quite nontrivial 
since the center of mass motion cannot be separated. 

Following Murthy et al. [5], we start with the iV-particle Hamiltonian (in 
the scaled variables) as given by 

^ ^ r2 ^ r--v. 2 

2i^ = -Ev? + ^iE^ + ^^ E ^^^-7= (33) 



where Xij is as given by eq. (p2D while gi and g2 are dimension-less coupling 
constants of the two-body and the three-body interactions respectively. 
It is easily shown that 



N I 

Mx^, y^) = E\ X,, 1^ exp{- 2 \Eo\J2 r') (34) 

i<j y i 

is the exact ground state of the system with the corresponding ground state 
energy being (in units oime^/h'^) 

^0 = -2[gN{N-l) + N-l]' ^^^^ 
provided gi and g2 are related to g by 

9i = gig - 1), g2 = g^- (36) 

It may be noted that this ipo is regular for g > which implies that gi > 
-1/4,^2 >0. 
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There is a neat way of proving that we have indeed obtained the ground 
state by using the method of operators [14]. In particular, let us define the 
operators 

N 
N 

Q,. = p.-.9i:|---||5 (37) 



and their Hermitian conjugates Q1^^ and Q^. It is easy to see that the Q's 
annihilate the ground state (0) i.e. Qx{ipo = Qyi'^o = 0. Further, the 
Hamiltonian (|33D can now be recasted in terms of these operators as 



+ Eo (38) 



where Eo{= —a'^/2) is as given by eq. (^). Since the operator on the right 
hand side is positive definite and annihilates the ground state wave function 



(|34D hence Eq must be the (bosonic) ground state energy of the system. Note 
that as in the other two- ( and multi-) dimensional problems, we are unable 
to obtain the fermionic ground state of the system analytically. In fact, we 
do not know of any A^-body problem in two or more space dimensions whose 
fermionic ground state has been analytically obtained. 

3.1 A class of Excited States 

We now show that a la the oscillator case, even in our case, a class of excited 
states can be obtained analytically. To that end, we consider the ansatz 

N I 



i'ixi^yi) = n I 1^ exp{-a J2^'i)^i^i^yi) (39) 



i<j 



On using eqs. (|33| ) and (|39D it is easily shown that ct^ = —2E while (j) satisfies 
the eigenvalue equation 



1 , ■ Vi A d d 



(40) 
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where 

A = a[gN{N -1) + N (41) 

The exact solutions of this differential equation are best studied in terms of 
the complex coordinates z = x + iy,z* = x — iy, and their partial derivatives 
9^ i = M- and d* ^ = + In terms of these 

coordinates the differential equation for take the form 



Ej ZjZ* ZjZ* i^j ZiZj ZjZi 



0. (42) 



It is worth remarking that (p is also an eigenstate of the total angular momen- 
tum operator i.e. L<l> = /$. We now classify some exact solutions according 
to their angular momentum, 
(a) Z = Solutions 

Let us define an auxiliary parameter t by t = 2ay/^~ZiZ* and let = 4>(t). 
On using this in eq. (^), the differential equation for reduces to 

t(f)"{t) + {b- t)(t)'{t) - a = (43) 

where h = 2N{N — l)g + 2N —1 , and a = A/a. The allowed solutions are [15] 
confluent hypergeometric functions, 0(t) = M{a,b,t). Thus the polynomial 
solutions are obtained when a = —rir with being positive integer. Here, 
the subscript r denotes radial excitations. The corresponding eigenvalues are 



2[nr + gN{N -1) + N 



Note that all these states have zero angular momentum, 
(b) / > Solutions 

Let tz = J2i and let = (pitz)- All the mixed derivative terms in eq. 
(^2p drop out and we obtain 

2atz<i)\tz) + A0 = (45) 

whose solution is (pitz) = where m = —A/2a. Note that 0(t^) is also an 
eigenfunction of the angular momentum operator with the eigenvalue / = 2m 
. Hence the energy eigenvalues in this case are given by 

" "2[(7iV(iV-l) + iV + /-i]2 ■ ^^^^ 
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(c) / < Solutions 

Similarly let tz = and = <P{tz*)- In this case (p satisfies the 

differential equation 

2atz*(f)'{tz*)+A(f)^0 (47) 

whose solution is (f){tz*) — f^* where m — —A/2a. In this case (f){tz*) is 
an eigenfunction of the angular momentum operator with the eigenvalue 
I — —2m < 0. Thus the energy eigenvalues in this case are given by 

" ~2[gN{N -l) + N ■ ^^^^ 

(d) Tower of Excited States 

We can now combine solutions in cases (b) and (c) (with nonzero I) with 
the solutions in (a) and obtain an even more general class of excited states. 
For example, let us define 

(/.(^„<) = 0i(t)(^2(t.) (49) 

where 0i is the solution with / = 0, while 02 is the solution with / > 0. The 
differential equation for 0i is again a confluent hypergeometric equation and 
the energy eigenvalues are given by 

" ~2K + ^iV(Ar-l) + iV + /-i]2 

One may repeat the same procedure to obtain exact sohitions for a tower of 
excited states with / < 0. Combining all these states, it is then clear that 
the exact energy eigenvalues may be written in the form 

^""^'"^ " ~ 2K + 9N{N - 1) + iV + 2 I m I -\Y ^^^^ 

We thus see that for both I > and Z < 0, one has a tower of excited 
states. Actually, the existence of a tower is a general result applicable to all 
excited states of which the exact solutions shown above form a subset. In 
particular, following the arguments given in Bhaduri et al. [6] , let us separate 



the coordinates (xj, yi) into one radial coordinate t = 2aJj2i ^1 ^ above and 
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2N — 1 angular coordinates collectively denoted by Vti say. Then eq. (^) 
can be expressed as 

t(t)"{t) + [2A^(A^ - 1)^ + 2A^ - 1 - t\(t)'{t) -[- + —T](f) = (52) 

where T = D2 + gDi. Here is an n'th order differential operator which 
only acts on functions of the angles On further factorizing 

(t>{x,,y,)=R{t)Y{n;) (53) 

where y is a generalized spherical harmonics on the (2N-l)-dimensional 
sphere 5"^^"^ it is easily seen that Y is an eigenfunction of the operator 
T. Let us say that the corresponding eigenvalue is A (of course it must be 
admitted here that the hard part of the problem is to find A). On further 
writing R{t) = t^R{t), it is easily shown that R{t) satisfies a confluent hy- 
pergeometric equation 

tR!'{t) + {b- t)R'{t) - akit) = (54) 

provided 

fi=l[N{N-l)g + N-lf + 4X] - [A^(A^ - 1)^ + - 1] . (55) 

Here b = 2N{N - l)g + 2N - 1 + 2fi while a = /i + A/ a. We thus get 
normalizable polynomial solutions of degree rir where a = —rir. Here rir = 
0, 1, 2... denote the number of radial nodes. Hence the energy eigenvalues are 
given by 

^ ~2K + <7iV(iV-l) + iV-i + /iP • ^^^^ 

It may be noted that fi is in general a complicated function of the coupling 
constant g. We thus see that for a given value of /x (which in general is 
unknown) there is an infinite tower of energy eigenvalues for which (—E)~^^'^ 
is separated by a spacing of one unit. This tower structure is a characteristic 
of the A^-body potential as given by eq. (|l]) or eq. (^) and will be present in 
any A^-body problem in Z)- dimensions. Note that a similar tower structure 
also exists in the case of any A^-body problem in D-dimensions if they are 
interacting via an oscillator potential except that in that case, for a given 
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value of /X, there is an infinite tower of energy eigenvalues separated by a 
spacing of two units. A la Bliaduri et al. [6], in our case too, the tower 
structure and the angular momentum are useful in organizing a numerical or 
analytical study of the energy spectrum. In particular, note that, the radial 
quantum number n^, and the angular momentum /, are integers, and hence 
they cannot change as the parameter g is varied continuously. 

Before ending this discussion it is worth pointing out that the class of 
excited states can also be obtained if we add the iV-body potential (|) to our 
A^-body potential as given by eq. In particular, on proceeding from eq. 
(|39D and repeating the steps, instead of eq. P3) one now has 



= . (57) 



A S 

+ , + 



EjZ.Z* ^EjZjZ* 

On following the steps as given above it is then easy to show that the exact 
energy eigenvalues are now given by 

^"-"^ = "2K + 7 + iP ^^^^ 

where 

J = ([gN{N -1) + N + 2\m\ -1]^ + 6^) . (59) 



The corresponding energy eigenstates are also easily written down. 

It is worth pointing out that if instead we add the A^-body potential (|^) 
to the oscillator potential, then on following the steps as given in [6], the 
exact energy eigenvalues can be shown to be 

= 2n, + 1 + 7 . (60) 



3.2 The Two-Body Problem: Complete Solution 

As in the oscillator case [6] , we now explicitly show that the two-body prob- 
lem is integrable and exactly solvable when the two bodies are interacting 
via the A^-body potential as given by eq. (|l]). As in [6] this is best done by 
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going over to the hyperspherical formalism in two dimensions [16,17]. It may 
be noted that the two-body problem is quite non-trivial here since the center 
of mass cannot be separated. 

The two body Hamiltonian is given by (see eq. (|33|)) 

H = -1(V? + V^) - , ^ + ^M±iil (61) 

where X = Xiy2 — X2yi- Note that the three-body term is obviously not 
there in the two-body problem ! The two body problem is best solved in the 
hyperspherical coordinates. To that end, let us parameterize the coordinates 
ri, r2 in terms of three angles and one length, {R, 9, 0, ip) as follows: 

Xi + iyi = R{cos 9 cos (j) — i sin 9 sin (j))exp{iip) 

X2 + iy2 = R{cos 9 sin (p + i sin 9 cos (f))exp{iip) ■ (62) 

For a fixed R, these coordinates define a sphere in four dimensions with 
radius R{R^ = rf + r|) and 9, cf) and if) are in the interval 

- 7r/4 < ^ < 7r/4 , -7r/2 < < 7r/2 , -tt < ^ < vr . (63) 

The important point to note is that X = Xiy2 — X2yi = -y- sin(26') depends 
only on R and 9 and is independent of and ^p. As a result the two integrals 
of motion are the angular momentum operator L 

d 

and the supersymmetry operator 

d d d 

Q = i[x2T^ + y2T: xi- y^—] = -i—. (65) 

oxi oyi 0x2 oy2 ocp 

Thus with SUSY, the two body problem as given by eq. ( [HT| ) is integrable 
with the four constants of motion being the Hamiltonian H, the angular part 
of H, L and Q. 

It is easy to check that the bosonic ground state has the quantum numbers 
/ and q of the angular momentum and SUSY operators equal to zero. In this 
context it is worth noting that the eigenstates of the SUSY operator Q are 
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neither symmetric nor asymmetric, unless the eigenvalue g = 0. After finding 
a simultaneous eigenstate of H, L and Q, we can separate it into symmetric 
(bosonic) and asymmetric (fermionic) parts. Note that these symmetric and 
anti-symmetric parts are eigenstates of but not of Q. 

The two-body Hamiltonian in terms of the hyperspherical coordinates is 
given by 



2 



d' 3d A' 2 



_dR^ ' RdR ' R\^ R^sm\29) ^^^^ 

where the operator is the Laplacian on the sphere and is given by 

2 _ d 2sin(2g) d 1 . 

~W cos(2^) 86^ co??{2eyd^^ ''d^^W 

As in the oscillator case, the interaction term i.e. 2gi/ R^ sin^(2^) is indepen- 
dent of <p and if). Hence the operators L and Q commute with the Hamiltonian 
( |66D (since they also commute with the noninteracting [gi = 0) Hamiltonian). 
Hence, for all gi, we label the states with the eigenvalues of L and Q. Each 
of these states is four-fold degenerate under parity L — > —L, Q ^ Q and the 
Hamiltonian is also invariant under parity. Hence the states labeled by the 
quantum numbers {l,q) have the same energy as {—l,q). The Hamiltonian 
is also invariant under the discrete transformation ri — ri and r2 r2. 
Note that under this transformation L ^ L and Q —Q. Thus the states 
labeled by the quantum numbers (/,g) have the same energy as (/, —q). In 
this way we obtain the four-fold degeneracy of the states. In fact we shall 
see below that the states with (/, q ) have the same energy as (g, /) since 
interchanging q and / leaves the differential equation invariant. As a result, 
one has in fact an eight-fold degeneracy for the levels for which | q \ and 
I / I are non-zero and different from each other. However, the degeneracy 
is only four-fold if | Z | = | g It^ 0. Similarly, there is a four- fold degeneracy 
between the states (±Z, 0) and (0, ±/) if / 7^ 0. It must be noted here that 
the degeneracy that we are talking about is a subset of the degeneracy of the 
noninteracting system. 

Let us now try to solve the eigenvalue equation Hip = Eip with H being 
as given by eq. (|66D. It is easily seen that if we write 

^(i?,^,0,V^) = F(i?)$(^,0,0) (68) 



16 



then the eigenvalue equation separates into angular and radial equations. In 
particular, the angular equation is given by 



[A' + ^%]'^ = /5(/5 + 2)$ (69) 
sm Zu 

where /?>—!, while the radial equation is given by 

F"{R) + ^F\R) + (2i? + I - ^^^)F(i?) = . (70) 
The radial equation is easily solved, yielding 

F{R) = R^exp{-J2 \ E \R)M{a, 6, 2J2 \ E \R) (71) 



where b = 2f3 + 3,a = 3/2 + f3 — l/^J2 \ E \ and M(a, b, x) is the confluent 
hypergeometric function. Demanding a = —rir, where is a positive integer, 
yields the bound state energy eigenvalues as 

^ = "2(n, + /3 + f)2- (^^^ 

The tower structure of the eigenvalues built on the radial excitations of the 
ground state is obvious from eq. ([721). 

It must be emphasized here that f3 is still unknown and has to be obtained 
by solving the angular equation ( |69D . We now note that the angular equation 
is in fact identical in the oscillator and our case and it has been analyzed in 
great detail by Bhaduri et al. 0. We can therefore borrow their results and 
draw conclusions about the value of (3 and hence the spectrum as given by 
eq. (0). Some of the conclusions in our case are 

1. On substituting 

^{9,(j),ip) = P{x)exp{iq(f))exp{il'(p) (73) 

where 

P{x) =1 X I" (1 - x)^(l + x)^e"'*'"(x) (74) 

with X = sin 9, b =\ I + q \ / 4, c = \ I — q \ / 4, a{a — 1) = gi = g{g — 1) 
and q being integers, it is easily shown that Q{x) satisfies the Heun's 
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equation 



;i - x^)e"ix) + 2[^-ib-c)-ia + b + c+l)x]Q'{x) 

+ [ ^ -(a + 6 + c+-)=^ + 2 ^ ^ ^ ]e(x)=0 (75) 



whose solutions are characterized by the so-called P-symbol |jT8|. Note 
that unlike the hypergeometric case, Heun's equation has four regular 
singular points. 

2. In case h = c (i.e. either / = or/and q = 0), the Heun's eq. (|75|) is 
exactly solvable [6] with corresponding P being 

/3 = 2m + 2a + 4b, m = 0,1,2, .... (76) 

Thus in this case E"^^"^ varies linearly with a (as is the case for the 
exact solutions of the many-body problem). As pointed out in [5], this 
is an example of a conditionally exactly solvable (CES) problem [ p!9[] . 



3. There are quasi-exactly solvable (QES) pO[ polynomial solutions of 



degree p{p > 1) in case there is a specific relationship between a, b and 
c . As a result, P and hence E~^^'^ is nonlinear in a. In particular, in 
this case /3 is given by 

/? = 2a + 26 + 2c + 2p . (77) 

It is worth emphasizing that the equation is exactly solvable at an 
infinite number of isolated points in the space of parameters (a, 6, c). 
These are isolated points because if we vary a slightly away from any 
one of them, the equation is not exactly solvable since b, c can only take 
discrete values and hence cannot be varied continuously. 

4. Bhaduri et al. [6] have done a detailed numerical, perturbative and 
large-^f analysis of /3 as a function of g both when Q < g < 1 and also 
for large g. Their analysis is also valid in our case. For example, their 
Figs. 1 and 2 are also valid in our case except that in our case, Fig.l 
represents a plot of /? + 2nj. + 2 as a function of g while Fig. 2 represents 
a plot of /? + 2 — 2(7 vs. g. Most of their discussion goes through in our 
case with this obvious modification. 
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The two body problem is also completely solvable in case we add the 
two-body potential as given by eq. to the Hamiltonian as given by 
eq. (pTD- In this case, the entire discussion about the angular part as 
given above is still valid. The only modification is in the radial equation 
as given by (|7^). In particular, instead of ([70|) we now have 

F"iR) + ^F'iR) + i2E+^- ^^^^§^)FiR) = . (78) 
whose solution is 



F{R) = R^exp{-J2 \ E \R)M{a, b, 2 J2 | E \R) (79) 



1 

2{nr + 7 + 



E = - o.^ , , 3,2 • (80) 



where 6 = 27 + 3, a = 7 + 3/2-1/^2 | E \ with 7 = ^(1 + /3)2 + 5^ - 1, 
while M(a, b, x) is the confiuent hypergeometric function. 

6. If instead if we add the two-body potential (^ to the corresponding 
oscillator problem, then it is easily seen that the angular part is again 
unchanged while the solution to the radial part is given by 

En^ = 2n, + 7 + 1 (81) 

F{R) = K'exp{-R^/2)M{a, 6, R^) (82) 
where a = 1 + (7 - E)/2, 6 = 7 + 2. 



4 Novel correlations in Arbitrary Dimensions 

In the last section we have obtained the exact bosonic ground state with 
Novel correlation as well as radial excitations over it in case the iV-particles 
in two dimensions also interact through the A^-body potential. The purpose 
of this section is to generalize these arguments to arbitrary number {D) of 
dimensions. In this context it is worth recalling that in case the N particles 
are interacting via the oscillator potential, such a generalization has recently 
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been carried out by Ghosh [7]. Following him, let us consider the following 
many-body Hamiltonian in D-dimensions 



1 ^ 1 n 

TT V72 ^ L ^1 ^i2i3-iD 



i=i y2^j=irj ^ R ^ hi2...iD 

2 ^ Piii2—iD-^iih—jD 

where it! denotes the sum over all the indices from 1 to with the restriction 
that no two indices can have the same value simultaneously. Here 

Thus the Hamiltonian has D-body and {2D — l)-body interactions only. 
Note that both Qi2i3...jo and Piii2...iD antisymmetric under the exchange 
of particle coordinates. Further, Pi^i2...iD vanishes when the relative angle 
between any two of the particles is zero or tt. 

For simplicity, let us first discuss the D — ?, case is some detail and then 
merely quote the results for the general case of D-dimcnsions. It is easily 
seen that the exact bosonic ground state in 3-dimensions is given by 



^0= X{\P^,k\'exp{-a Y.^}) (85) 

R V i 



where Eq = provided 



5i = |(f-l), 52 = ^. (86) 
The corresponding ground state energy Eq is given by 

" ~2[\gN{N-l){N-2) + ^-f-\]^ ' ^^^^ 

The fact that this is indeed the ground state is easily shown as follows. 
Consider the following annihilation operators in three dimensions (note that 
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n = Xii + yij + Zik) 



A. 



^ s 



:q 



jkjxi 



Pz 



2 ^ 



Pijk 




{^jk)yi 


iaui 


Pijk 




(Qjr'fc)zi 


iazi 


Pijk 





(88) 



where S denotes the sum over all the repeated indices from 1 to with the 
constraint that any two indices cannot have the same value simultaneously. 
It is easily shown that the Hamiltonian (|8^) in 3-dimensions can be written 
down as 

1 ^ / \ 

E UtA^. + + AtA^. + (89) 



H 



i=l 



Further, it is easily checked that the annihilation operators A's annihilate 
the wave function ipQ as given by eq. (|^) and hence ipQ is indeed the ground 
state with the ground state energy being Eq as given by eq. (|57D. 

The excited states can also be obtained in our case [7]. To that end we 
write 

i}){xi, Ui, Zi) = iJoixi, yi, Zi)(t){xi, yi, Zi) (90) 

where ipQ is as given by eq. (|85D except that a is now given by = —o? 12. 
On substituting this ansatz in the Schrodinger equation for the Hamiltonian 
13D, it is easily shown satisfies the equation 



1 ■tr^ —o 9 sr^ Qjk ■ Vi 



2 E 2 E 



P 



a 



ijk 



E 



ri ■ Vi 



A 







(91) 



where 



A 



a 



|iV(iV-l)(iV-2) + ^-i 



On assuming that is a function of t = ^ayj^^h it is easily seen that eq. 
ITF) for (j)(t) reduces to the confluent hypergeometric equation 



(92) 



t(f)"{t) + {b-t)(j)'{t) -a(j){t) = 



(93) 
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when a = A/a^ b = 2{A + l)/a. Hence the admissible normahzable solutions 
of eq. (|93D are M(a = —nr,b,t) with the corresponding energy eigenvalues 
being 

^ "2K + fiV(iV-l)(iV-2) + f • ^^^^ 

Note again the tower structure and the fact that {—En^Y^"^ is linear in g for 
the exact solutions. 

We can also obtain the exact solutions in case we add the A^-body po- 
tential to the Hamiltonian (PB]). In particular the energy eigenvalues are 
now given by 

= -2ktW 

where 

,= (l^E±lEiIl^ME^.^^^sf\ (96) 

The corresponding eigenfunctions can be easily written down, if instead one 
adds the A^-body potential (H) to the oscillator potential, then the energy 
eigenvalues are given by 

= 2n, + 1 + 7 . (97) 

The corresponding eigenfunctions can again be easily written down. 

The generalization of the above results to D-dimensions is straight for- 
ward. In particular, it is easily shown that the exact bosonic ground state and 
the tower of excitations over it in the case of the Z)-dimensional many-body 
Hamiltonian (|83D are given by 



V'n = n I Piih-iD \' M{a = -rir, b, t)exp{-a r?) (98) 

R V i 



where 



provided 



2 2[nr + !^ + gD{^CD)~\Y ^^^^ 



^-((z^)'^-(z^[(z^-^] « 
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so that 

9=^^^[l±^l+Ag,] (101) 

It is interesting to note that whereas the relation between gi and g2 and also 
between g and gi is Z)-dependent, the allowed ranges of gi are independent 
of D i.e. in any dimensions, whereas the solutions in the upper branch are 
regular for gi > — | , in the lower branch the regular solutions are only allowed 

in the limited range — | < (71 < 0. Here a = A/ a^h = ^^^"^-^ where 

A = a(|[iV + 2^(^Cz,)]-^)-l. (102) 

The exact solutions are also obtained if we add the A^-body potential (|^) 
to the Hamiltonian (^) and the exact energy eigenvalues are given by 

where 

l={[^ + gD{^Cn)-lf + 5') . (104) 

Similarly, if we add the A^-body potential (^ to the Hamiltonian (|83|) but 
with the oscillator potential then the exact solutions are given by 

En, = 2nr + l + -i . (105) 

The corresponding eigenf unctions can be easily written down in both the 
cases. 



5 Calogero-Sutherland Type Models in Higher 
Dimensions with N body Interaction 

In Sec. II, we have considered one possible generalization of the Calogero- 
Sutherland type models in higher dimensions. The key point there was to 
have a long-ranged three-body interaction term. This is over and above 
the long ranged two-body interaction term which is present in one dimen- 
sion. Only then it was possible to obtain a class of exact solutions includ- 
ing the bosonic ground state. Another possible generalization of Calogero- 
Sutherland model to higher dimensions was considered recently by Ghosh [8] 
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when he introduced two models with purely two-body long ranged interac- 
tions and in both the cases he was able to obtain the exact bosonic ground 
state and radial excitations over it in case the iV-particle also interact via a 
(one-body) oscillator potential. In this section we shall show that the exact 
ground state and radial excitations over it can also be obtained in case the 
oscillator potential is replaced by the A^-body potential as given by eq. (|l|). 
Following Ghosh [8], let us consider the Hamiltonian 

^ = -o E - + v,ip) + v^iP) + WM (106) 



2 k=i JTk rf 



where 





r 


'k 


|2(/3-l) 


(|r. 






|/3)2 







(/3-2) 


(|r. 


|/3 




1^) 



^ k+i 
g2 I |2(/3-l) 

wm - y«to-^.|:, (|,|.-|^;iVi^-k>H 

where g is a dimension-less constant. Note that for (3 = 1 and 2 the three 
body interaction term vanishes. Thus even though we give results for 
arbitrary positive values of /5, it must be noted that one has long ranged 
two-body interaction alone only if /3 = 1,2. 

We first note that the Hamiltonian ( |106| ) can be written as 

H = WK■^^ + Eo (108) 
^ i 

where the annihilation operators Aj are given by 

I V' \^ ^ z(yT ' 

ijLj \\ i \ I i I j 

while the ground state energy Eq{= — a^/2) is 

^° = "2[i^ + fiV(iV-l)-ip- 
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The fact that Eo is the bosonic ground state energy is verified by noting that 
each of the operators Aj annihilate the ground state wave function 



^0 = n(l 1^ - I fyexpi- 2 \ E, \ Y^r}) . (Ill) 

It is straight forward to obtain the radial excitation spectrum. In partic- 
ular, it is easily shown that the corresponding exact eigenf unctions are 



^ = lid 1^ - I r, |^)SM(a = ~nr,b,t)exp{- 2 \E\Y.rf). (112) 

i<j y i 

with the corresponding eigenvalues being 

^ "2K + f^ + fiv(iv-i)-i]2 • ^^^^^ 



Here t = 2a^Jj2ir],a = A/a,b = 2{A + l)/a while E = -a^/2 and A = 
a[ND/2 + gpN{N - l)/2 - 1]. 

One can also obtain the spectrum in case one adds the A^-body poten- 
tial (^) to the Hamiltonian ( |106| ). In particular, it is easily shown that the 
corresponding energy eigenvalues are 

where 

The corresponding eigenfunctions can also be easily written down. Similarly, 
in the oscillator case too, the spectrum can be written down in the presence 
of the A^-body potential (^. 



6 Complete Bosonic spectrum of an A^-body 
problem in D dimensions 

Over the years, the exact solution of an A^-body problem in three (and even 
arbitrary) space dimensions has attracted considerable attention because of 
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its obvious implications in several areas. So far, the only A^-body problem 
which is completely solvable in two and higher dimensions is when A^-bosons 
are interacting pair-wise via harmonic interaction. The purpose of this sec- 
tion is to show that there is one more A^-body problem in two and higher 
dimensions for which the complete discrete spectrum can also be obtained. 
In particular I show that the full bound state spectrum can also be obtained 
in case A^-bosons are interacting via the A^-body potential as given by eq. 

(i- 

In particular, consider the Hamiltonian 

1 ^ 1 

H = -^J:^'.- ;^ . (116) 



i=l 



After the separation of the center of mass (which in this case moves as a 
free particle), the relative problem is usually discussed in terms of the Jacobi 
coordinates = 1,2, N — 1) defined by 

= (-^)^/^[j(r. + ... + r,)-r,^.i]. (117) 

However, it is more convenient to consider the problem in hyperspherical 
coordinates {(, u) of the (A^ — l)-dimensional vector ^. Here the hyper-radius 
C(0 < C ^ oo) is given by 

N-l 

= E c| (118) 

while is a set of [(A^ — 1)-D-1] angular coordinates. With this choice, the 
Hamiltonian takes the form 



2 



92 ^{N-1)D-1 d L\iu) 



^ (119) 



A^C 



where L'^{uj) is the "grand angular" operator whose eigen functions are the 
orthonormal hyper-spherical harmonics Y[l]{uj) with eigenvalues —k[k + {N — 
1)D — 2]. Here k = 0,1,2,... is the grand angular quantum number. It is thus 
clear that if we write the eigenfunction ip in the form 

^(C,a;) = 0(C)y[L](^) (120) 
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then in the ( variable one has essentially a one dimensional Schrodinger 
equation 



d'^ ^{N -1)D ~1 d k[k+ {N -1)D -2] ^ 



(121) 

for the Kepler problem in {N — l)Z}-dimensions. The corresponding discrete 
energy energy values and eigen functions are given by 

1 

^ "2iVKTI+^^^i^ ^^^^^ 

m = y'exp{-y/2)L'^(y) (123) 



where a = {N - l)D/2 + k - l,y = 2^2 | E |C and = 0, 1, 2, ... is the 
radial quantum number. 

It is worth pointing out that in eq. (|119|) , apart from the iV-body potential 
the only other potential for which the Schrodinger equation can be solved 
for all values of k is the two-body harmonic oscillator potential. Besides, in 
both the cases one can also add the A^-body potential and the problem 
is still analytically solvable. In particular, instead of ( |122| ) the spectrum is 
now given by 

1 



2iVK + i + A]2 ^^^^^ 



where 



A = Va2 + 52 . (125) 

The corresponding eigenfunctions can also be easily written down. Note 
however that the degeneracy in the spectrum is now very much reduced. 

Our A^-body potential is in a way richer than the oscillator potential 
in that here one not only has an infinite number of negative energy bound 
states, but one also has positive energy scattering states. In particular, for 
£^ > 0, the solution to the Schrodinger eq. ( |121| ) is given by 

UQ = exp{ipC)C'F(k + ~ ~ ^ - -, (iV - l)D + 2K-1, -2^k) 

^ (126) 
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where E = p'^/2. Hence the phase shifts for this problem are 



p2i<5p ^ _^ 1 Ei (127) 



Unfortunately, we are unable to draw any simple conclusion about the A^- 
particle scattering from this expression of the phase shift. 



7 Summary and Open Problems 

In this paper, we have discussed several A^-body problems in two and higher 
dimensions and have provided support to the conjecture that whenever an A^- 
body problem is (partially) solvable in case the iV-bodies are interacting by 
the harmonic forces, then the same problem will also be (partially) solvable 
in case they are interacting by an A^-body potential as given by eq. ([^) (or 
(^). Our conjecture is based on the following simple observation. In all the 
many-body problems with harmonic forces discussed in [2-8,13,17], after the 
short distance correlation etc. is taken out, the problem essentially reduces to 
that of harmonic oscillator in dimensions higher than one. Now the only two 
problems for which the bound state spectrum can be analytically obtained 
for all partial waves in two and higher dimensions are the oscillator and the 
Coulomb problems. This strongly suggests that those A^-body problems for 
which the Schrodinger equation essentially reduces to Coulomb problem in 
two and higher dimensions, after the short distance correlation part etc. is 
taken out, must also be (partially) solvable. While we have no proof for 
this conjecture, we have not come across any counter example either. We 
have also studied other models in one and two dimensions with this A^-body 
potential [9,10] which also provide support to this conjecture. It is clearly 
necessary to examine other solvable many-body problems with the harmonic 
forces and check if our conjecture is true in those cases or not. 

Another well known fact in non-relativistic quantum mechanics is that 
both the Coulomb and the oscillator problems in two and higher dimensions 
are also exactly solvable for all partial waves if one adds a potential of the 
form ^^5^/2mr^ to either one of them. Based on this observation, we have 
conjectured that for the A^-body problems with the oscillator or our A^-body 
potential as given by eq. (P (or (|)), one can also add the A^-body potential 
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as given by eq. (|D (or by (||)) and the problem is still (partially) solvable 
but now the degeneracy in the spectrum is much reduced. While we have 
no proof for this conjecture, we as well as Gurappa et al [13] have provided 
support to this conjecture through several examples. Besides, we have not 
come across any counter example either. 

Apart from the oscillator and the A^-body potential (|lD, are there other 
potentials for which any of these A^-body problems in higher dimensions can 
be (partially) solved? While we have no definite answer, it appears that 



all other potentials will only be quasi-exactly solvable since we know of 
no other potentials in non-relativistic quantum mechanics in two and higher 
dimensions for which the bound state spectra can be analytically obtained 
for all partial waves. 

There are several common features in all the A^-body problems that we 
have studied in this paper and before [9,12] as well as those studied with the 
oscillator potential. Some of these are as follows. 

1. By now several A^-body problems exist in two and higher dimensions 
for which a class of eigenstates including the bosonic ground state have 
been exactly obtained in case they also interact either via the oscillator 
or the A^-body potential (||). All these problems have long ranged two- 
body and in most cases also the long ranged three-body interactions. In 
the special case of one dimension, there are only long ranged two-body 
interactions and the complete bound state spectrum has been obtained 
in both the cases [2,10]. 

2. In all the A^-body problems with the oscillator or the A^-body potential 
studied so far, one always obtains the tower of states characteristic 
of the oscillator or the Coulomb problem as the case may be. We 
conjecture that this will be a general feature of any many-body problem 
with either of these forces. While we have no general proof, we have 
not come across any counter example either. 

3. All the A^-body problems in two and higher dimensions (that have been 
studied so far) with either the oscillator or the A^-body potential are 
only partially solvable in the sense that while the bosonic ground state 
and radial excitations over it are known, the full bosonic spectrum is 
still unknown (except the one studied in the last section or its oscillator 
analogue). It is worth pointing out that the same is also true in the 
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case of iV-anyons interacting either via the oscillator or via our A^-body 
potential [17]. Further, for none of these problems, the ground state 
of A^-fermions is analytically known. I conjecture that the same will 
be true for any other A^-body problem in two and higher dimensions 
interacting via either of these potentials. 

4. If anyon example is any guide, then it would seem that in the case of 
any A^-body problem in two and higher dimensions with the oscillator 
potential, a part of the spectrum will be linear in the coupling constant 
(say Ad in Sec. II) while a part of the spectrum will be nonlinear in 
the coupling constant. On the other hand, in the case of our A^-body 
potential (||), a part of (— will be linear and a part will be 
nonlinear in the coupling constant. So far, in all the examples discussed 
in the literature, one has only been able to obtain the linear part of E 
(or (— E")^/^) in the case of the oscillator (or A^-body) potential. While 
it is not clear if one has obtained the full linear spectrum in all the 
cases or not, it is certainly true that so far exact solution has not been 
obtained for even one nonlinear state in any of the A^-body problems 
in two and higher dimensions with either of the potentials. I believe 
that if one can obtain exact solution for even one nonlinear state in any 
of these problems, it will be a major breakthrough. 

It must be added here though that the oscillator potential is perhaps more 
useful than our A^-body potential since whereas the clustering property holds 
in the case of the oscillator potential it does not hold in the case of the A^-body 
potential. On the other hand, the A^-body potential is in a way richer than 
the oscillator potential in that unlike it, one has both bound and continuous 
spectrum in this case. The most important and of course difficult problem is 
to find some physical application of the A^-body potential. Hopefully, some 
application will be found in the near future. 
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